Multiband model for penetration depth in MgB 2 
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The results of first principles calculations of the electronic structure and the electron-phonon 
interaction in MgB2 are used to study theoretically the temperature dependence and anisotropy of 
the magnetic field penetration depth. The effects of impurity scattering are essential for a proper 
description of the experimental results. We compare our results with experimental data and we 
argue that the two-band model describes the data rather well. 
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The electronic structure of the recently discovered 
superconductorty MgB2 is now rather well understood 
and the superconductivity may be asa^ed to the con- 
ventional electron-phonon mechanism.tmcra The Fermi 
surface consists of two three-dimensional sheets, from 
the 7r bonding and antibonding bands, and two nearly 
cylindrical sheets from the two-dimensional a bands. 
The qualitative difference between the 2-D a- and the 
3-D 7r-bands in connection with the large disparity of 
the electron-phonon interaction (EPI) for the different 
Fermi surface sheets_sjjftefis±ed a multiband description 
of supercondjictivityEu'Suilj Recent reports on quantum 
oscillationstil provided not only important information on 
the electronic structure near the Fermi level but it also 
probed directly the disparity of the EPI in the n- and 
tr-band systems. The excellent agreement of the calcu- 
lated EPI with the de Haas - van Alphen mass renormal- 
ization cleacly confirms the basic assumption of the two- 
gap model.!! 2 ] Further experimental support of this, mod el, 
comes from STM and point-contact spectroscoayjl3llil't3 
high-resolutioa, photo-emission spectroscojawij Raman 
spectroscopy ,113 specific heat measurements Mp and stud- 
ies of the magnetic penetration depth. Ej : E30 

There is still some debate concerning the applicabil- 
ity of the multiband description to MgB2, in particular 
since some tunneling measurementsE3 show only a sin- 
gle gap with a magnitude smaller than the BCS value of 
A = 1.76 T~. A recently proposed multiband scenario 
for tunnelingB in MgB2 explains the reason for the differ- 
ences in the observed tunneling spectra and thus helps 
to settle this debate. 

A similar discussion has been going on concerning the 
penetration depth. The measured magnetic penetration 
depth shows a large variety of behavior (see Table |). 
In order to interprete these results, a microscopic model 
is reqaked. In this paper we shall use the multiband 
modeE'B'B to calculate the temperature dependence and 
the anisotropy of the penetration depth using the Eliash- 
berg formalism and the results of first principles elec- 
tronic structure calculations. 

Generalization of the BCS theory to the multiband 



model was first suggested in Refs. 



and 



and 



inNb doped SrTi0 3 .B 
!J suggested a two-band 



has been observed experimental!; 
More recently, Kresin and Wolf 
model in the strong-coupling regime for describing the 
properties of high T c superconductors. Strong-coupling 
two-band-model calculations of the microwave response, 
and in particular the penetration depth, were performed 
in Refs. ^8| and |39| . 

The penetration depth of the magnetic field \l.oi/3 in 
the local (London) limit is related to the imaginary part 
of the optical conductivity by 
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where a, (3 denote Cartesian coordinates and c is the 
velocity of light. If we neglect strong-coupling effects 
(or, more generally, Fermi-liquid effects) then for a clean 
uniform superconductor at T — we have the relation 
X L , a p = c/ujf, where {ujf) 2 = 8ne 2 £ kj S(e^)v F v F is 
the squared plasma frequency and v F the a-component 
of the Fermi velocity. Impurities and interaction effects 
drastically enhance the penetration depth, and it is there- 
fore suitable to introduce a so called 'superfluid plasma 
frequency' oj s J a p by the relation u^ a/3 = c/X L , al3 . It has 
often been mentioned that this function corresponds to 
the charge density of the superfluid condensate, but we 
would like to point out that this is only the case for non- 
interacting clean systems at T = 0. 

In the two-band model we have J-hi^standard expres- 
sion (neglecting vertex corrections)' 38 ' 
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where Co{n) = u} n Z(oj n ) and A(u) n ) = A4u n )Z(u} n ) are 
the solutions of the Eliashberg equationsBEl and the cal- 
culated plasma frequencies for the a and tt bands are 
given in Ref. [| The Eliashberg equations were solved 
numerically as described in Refs. ^ and |^. The influence 
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TABLE I: Penetration depth measurements by different methods and groups (MW=microwave, /iSR=muon spin relaxation, 
RF=radio frequency, FIR=far infrared optical spectroscopy). Values for the estimated London penetration depth Al(0), 
superfluid plasma frequency oj^ , temperature dependence AA(T), superconducting gap values Ao and ratios 2Ao/IcbT c are 
shown. 
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AtSR+AC 
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2.33 
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RF 


160±20 


1.24 


BCS 
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FIR 


300 
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2.5 
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MW (T C =37.9K) 


100 


1.98 












MW (T C =26K) 


1200 


0.165 
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MW (T c =39K) 


110 


1.8 


BCS, T > 5K 


3.8 


2.26 






MW (T c =36K) 


115 


1.72 


BCS, T > 5K 


3.2 


2.06 







of impurities is incorporated into the model by including 
shifts of the gap function A°(w„) and the renormalization 
factor Zf(ui n ) 

A 4 - A° + £ 7u A°/2^ + (A°) 2 , 

3 

Zi{u n ) - ZPK)+X!7072^ + (A?) 2 
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in the Eliashberg equations. Intraband scattering does 
not change T c and the gap values (Anderson's theorem), 
but influences strongly the penetration depth. 

Before we start discussing the exact solutions to the 
Eliashberg equations we present a simplified model con- 
sisting of two independent BCS superconducting bands 
with different plasma frequencies and different gaps (and 
consequently different T c 's). In spite of the fact that this 
model is clearly an oversimplification, it captures quali- 
tatively most of the observed behavior. In this model the 
band with the larger T c has the smaller plasma frequency 
(see Fig. |). 

For a clean system the resulting inverse squared pen- 
etration depth is the sum of the 'superfluid plasma fre- 
quencies' (solid line). The kink in Fig. [j] is an artefact of 
the simplified model which will be smoothed out by inter- 
band scattering, but an inflection point in the tempera- 
ture dependence of the penetration depth remains which 
has been observed experimentally. The low temperature 
dependence is determined by the band with the small- 
est gap, whereas the high temperature behavior results 
from the band with the larger gap. This is in accordance 
with the temperature dependence for low temperatures 
observed in some experiments. If the superconducting 
band with the smaller gap will be 'overdamped' due to 
impurities, then the penetration depth is only determined 
by the other band and it will show a BCS temperature 




FIG. 1: Temperature dependence of the penetration depth 
in the model of two independent BCS superconducting bands 
(dashed and dotted line) with different superconducting gaps. 
The resulting penetration depth (solid line) clearly shows a 
non-BCS temperature behavior. The low temperature behav- 
ior will be dominated by the band with the smaller supercon- 
ducting gap. 



dependence, which has also been observed in some ex- 
periments. 

For a proper understanding of the observed physical 
behavior of MgB2 it is important that impurities are 
taken into account properly. Recently, the influence of 
impurities_on the two-gap superconductivity has been 
discussed.Ej Using the same arguments we shall discuss 
two cases: (a) The clean case with scattering rates-^j/o- = 
7 7r =2meV as realized in low-resistivity dense wirescll and 
(b) the dirty case with 7 CT =54meV and 7 7r =1.2eV. The 
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FIG. 2: The calculated temperature dependence of the pen- 
etration depth for the clean and dirty case (as defined in the 
text) in the a&-plane and along the c-axis, as well as a BCS 
curve corresponding to a single band case. For comparison ex- 
perimental data from microwave experiments on single crys- 
talsr-fA),LJ oriented films (o)E3 and fiSR data on polycrystals 
(■)li£l are shown too. 
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FIG. 3: The anisotropy of zero-temperature penetration 
depth vs impurity scattering rate in the 7r band. The in- 
set shows the temperature dependence of the anisotropy for 
the clean limit. 



values for the scattering rates in the dirty cassias well 
as negligibly small interband scattering rates,e2l are-, in 
accordance with the results on high resistivity films." 

Exact calculations, i.e. solving the Eliashberg 
equations for the effective two-band model with pa- 
rameters derived from first-principle electronic struc- 
ture calculations, have been carried out for the clean 
and dirty cases for a magnetic field along the c- 
direction or in the a6-plane. The results obtained can 
be presented in the form of the effective superfluid 
plasma frequency, wif . Fig. || displays the calculated 

(oj s /(T)/oj s /(T - 0)) 2 = (A L (T -» 0)/X L (T)) 2 as a func- 
tion of reduced temperature. First, we shall discuss 
the temperature dependence of A£ , when the magnetic 
field is oriented exactly along the c-axis (this means that 
screening currents run in the afr-plane). In the clean 
case the situation is similar to the model discussed above 
(Fig. |l|). \f^(T)^ 2 has an inflection point and the low 
temperature behavior is determined by the band with 
the small gap A T . In the dirty case the conductivity 
in the 7r-band is strongly suppressed. This means that 
the screening currents in the afr-plane are determined by 
cr-band with a BCS-like temperature dependence with 
a large gap A a . For the intermediate case the tempera- 
ture dependence of \°£(T) is between these limiting cases. 
One can even have situations with a nearly linear depen- 
dence in some temperature interval, as may be seen in 
Fig. ||. Experimepial data from microwave experiments 
on single crystalsE£l .and oriented fllmslS as well as /iSR 
data on polycrystalsta are shown for comparison in Fig. ||. 

As may be seen from Eq. (0), the penetration depth 
in the c-direction in the clean case is only determined 
by the 7r-bands because of the very small plasma fre- 
quency of the cr-band in this direction. It is interesting 



to notice that the temperature dependence of the pen- 
etration depth in the dirty c-axis case follows the clean 
ab-direction case. This is because both the tr-channel is 
blocked by the small plasma frequency and the 7r-channel 
is blocked by impurity scattering. However, the abso- 
lute values of the penetration depths differ by a factor of 
about 8 in these two cases. One may see that an inflec- 
tion point is present in the temperature dependence even 
for this single-band contribution, because of the induced 
superconductivity at higher temperatures in the 7r-band. 

The corresponding London penetration depths at T=4 
K have the values \fffi=39.2 nm, Af e c a ™=39.7 nm, 

^hlb ~ 105.7 nm, and A^ iy =316.5 nm. One may ob- 
serve that even in the clean case the value of w s J a p (T — > 
0) ~ 5eV from Eq. || differs from the total plasma fre- 
quency ~ 7 eV as a consequence of strong-coupling effects 
due to electron-phonon interaction. 

Table |l] summarizes the experimental information on 
the penetration depth Xl(T — > 0), the temperature de- 
pendence of the penetration depth AA l (T) and the esti- 
mated superconducting gap obtained by different exper- 
imental methods and groups. Our theoretical values of 
the penetration depths for the clean case are smaller than 
the smallest experimental value. On the other hand the 
values for the dirty case are in reasonable agreement with 
experiment. Nearly all measured penetration depths fall 
within our limiting cases (clean and dirty) , and especially 
the BCS-like behavior at lower temperatures observed 
in experiment, reflecting the 7r-band contribution, is in 
agreement with our theoretical calculations. 

It is well known that in a clean superconductor the 
low temperature penetration depth is independent of the 
superconducting gap. In this case, the anisotropy is only 
determined by the ratio of the plasma frequencies in the 
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ab-plane and c-directions. Hence, if the 7r-band is very 
clean the penetration depth is nearly isotropic due to the 
small difference in the effective plasma frequencies. This 
is shown in Fig. [| where one may also see that impu- 
rities in the 7r-band drastically enhance the anisotropy. 
In the inset of Fig. || the temperature dependence of the 
anisotropy is shown for the clean case. The reason for 
the strong variation with temperature is that for high 
temperatures the difference of gaps also contributes. A 
similar observation has recently been made in a weak 
coupling model.E3 

A final remark concerns the orientation of the magnetic 
field. According to estimates in Ref. || the u-band does 
not contribute to the electronic transport for angles with 
the c-axis larger than of the order of 1°. This implies that 
for larger angles the effective penetration depth is deter- 
mined by the 7r-band only. Only for angles approaching 
zero, the <r-band contributes and the penetration depth 
decreases towards the minimal value corresponding to the 
screening current flowing in the ab plane, namely 39 nm 
and 106 nm for the clean and dirty case respectively. 

The above mentioned considerations must be taken 
into account when interpreting the experimental data. In 
polycrystalline samples the penetration depth is mostly 
determined by the 7r-band and therefore practically 
isotropic and it is similar to the c-axis penetration depth 
in Fig. 0. Our calculations for the dirty case describes 



qualitatively well the data in Ref. |l9| On the other 
hand, the data for single crystals and oriented films cor- 
respond to our calculations of the a6-plane penetration 
depth, provided the magnetic field is oriented with an 
accuracy better than 1° along the c-axis. Therefore the 
data from Ref. ^6] and ||| are described by our calcu- 
lation for the clean case. Quantitative deviations can 
be attributed to a different impurity content and a pos- 
sible admixture of c-axis contribution. The temperature 
dependence of the specific electrical resistivity is not pro- 
vided in these papers however, which would be needed in 
order to estimate the impurity scattering rates. 

In conclusion, we have used the results of first princi- 
ples calculations of the electronic structure and electron- 
phonon interaction in MgE$2 to calculate the magnetic 
field penetration depth. The measured temperature de- 
pendence of the penetration depth is qualitatively well 
reproduced in a two band model with the same, set of 
parameters which was used to fit DC resistivity .El-El We 
predict strong dependence of the anisotropy of the pen- 
etration depth on impurity scattering. This anisotropy 
increases with increasing temperature. 
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